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Abstract We provide a phenomenological analysis of present experimental 
searches for local parity violation manifested through the Chiral Magnetic Ef- 
fect. We introduce and discuss the relevant correlation functions used for the 
measurements. Our analysis of the available data from both RHIC and LHC 
shows that the present experimental evidence for the Chiral Magnetic Effect 
is rather ambiguous. We further discuss in some detail various background 
contributions due to conventional physics, which need to be understood quan- 
titatively in order to draw a definitive conclusion about the existence of local 
parity violation in heavy ion collisions. 
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The theoretical study of topological solitons in field theories has a long his- 
tory. Quite generically, these objects arise as solutions to the classical equa- 
tions of motion for field theories due to the nonlinearity of the equations as 
well as due to specific boundary conditions. They are found in field theo- 
ries of various dimensions (2D kinks, 3D monopolcs, 4D instantons), and are 
known to be particularly important in the non-perturbativc domain where 
the theories are strongly coupled. For a recent review, see e.g. [Tj. 

Topological objects in Quantum Chromodynamics (QCD) are known to 
play important roles in many fundamental aspects of QCD [T|. For example, 
instantons are responsible for various properties of the QCD vacuum, such 
as spontaneous breaking of chiral symmetry and the Ua(1) anomaly (see 
e.g. [2j |3]). Magnetic monopoles, on the other hand, are speculated to be 
present in the QCD vacuum in a Bose-condensed form which then enforce 
the color confinement, known as the dual superconductor model for QCD 
confinement, which is strongly supported by lattice QCD calculations (see 
e.g. [H[S]). Alternatively vortices are believed to describe the chromo-electric 
flux configuration (i.e. flux tube) between a quark-anti-quark pair in the 
QCD vacuum which in turn gives rise to the confining linear potential (see 
e.g. reviews in [5JIE]). Some of these objects, such as monopoles [7] and flux 
tubes [5] , may also be important degrees of freedom in the hot and deconfincd 
QCD matter close to the transition temperature T c , and may be responsible 
for the observed properties of the so called strongly coupled quark-gluon 
plasma 0QI5]. 

Since the existence of such topological objects is theoretically well mo- 
tivated and their effects on the dynamics are deemed to be important, a 
direct experimental detection of such objects or at least of certain unique 
imprints by them, would be a highly desirable goal. This review will dis- 
cuss recent efforts and progress toward that goal, specifically in the context 
of rclativistic heavy ion collisions through the measurement and analysis of 
charge-dependent correlations. 



1.1 The Chiral Magnetic Effect in brief 

An interesting suggestion by Kharzeev and collaborators [HJ [12j UHl HH 
[T51 HB] on the direct manifestation of effects from topological objects is the 
possible occurrence of V- and CP-odd (local) domains due to the so-called 
sphaleron or anti-sphaleron transitions in the hot dense QCD matter created 
in relativistic heavy ion collisions. Imagine that in a single event created in a 
heavy ion collision the gauge field configurations in the space-time zone of the 
created hot dense matter experience a single sphaleron transition. As a re- 
sult this local zone acquires a non-zero topological charge which is parity-odd. 
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This non-zero topological charge, when coupled with light quarks through the 
triangle anomaly, induces a non-zero chirality for the quarks. In other words 
it generates an imbalance between left- and right-handed quark numbers, or 
a non-zero axial charge density. To be precise, there is no violation of parity 
at the interaction level, but rather a local creation of matter with non-zero 
axial charge density, which is a V- and CP-odd quantity. 

A concrete proposal for experimental detection is the so-called Chiral Mag- 
netic Effect (CME) [13]. The effect itself states that in the presence of external 
electromagnetic (EM) magnetic field B, a nonzero axial charge density will 
lead to an EM electric current along the direction of the magnetic field B: 

iv = ^ B (!) 

where fiA is the axial chemical potential associated with the non-zero axial 
charge density present in the system, and N c is the number of colors. This 
elegant relation is theoretically well established in both the weakly-coupled 
and the strongly-coupled regimes of the theory as will be discussed in several 
contributions to this volume. 

At first sight, it might seem that the above relation is violating parity: 
under spatial rotation and inversion the EM electric current jy transforms 
like a vector, while the magnetic field, B, transforms like an axial- or pseudo- 
vector. Therefore, the factor in Eq. (fT]) relating the two will have to be parity- 
odd. This is indeed the case, since /i^ that enters the above relation is a 
pseudo-scalar quantity which changes sign under parity transformation. Thus 
the CME relation, Eq. (TTJ), is invariant under parity transformation. However, 
in a region with nonzero, either positive or negative, (ia certain parity-odd 
obscrvables, e.g. the pesudoscalar quantity < jy ■ B >, may acquire nonzero 
expectation values. It is only in this sense that one may refer to it as "local 
parity violation" . 

In addition there is a complimentary relation, as one might have guessed 
from the "duality" by interchanging the roles of V (vector) and A (axial), 
that has been called the Chiral Separation Effect (CSE). The CSE refers to 
the separation of chiral (or axial) charge along the axis of the external EM 
magnetic field at finite density of the vector charge, for example at finite 
baryon number density |17[ 118] . The resulting axial current is given by 

Ja = ^#MvB (2) 

with the fj,y here being the baryon number chemical potential. Furthermore 
the combination of the two effects, CME and CSE, gives rise to an interest- 
ing propagating collective mode: the vector density induces an axial current 
which transports and creates a locally nonzero axial charge density, which 
in turn leads to a vector current that further transports and creates a lo- 
cally nonzero vector density, and so on. This is called Chiral Magnetic Wave 
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(CMW) [15] , just like Maxwell's electromagnetic waves represent the coupled 
evolution of the electric and magnetic fields. The CMW is a general concept 
that includes both the CME and CSE effects. It is robust in the sense that it 
takes the form of a collective excitation like the sound wave without relying 
on a quasi-particle picture. 

We end the general introduction with two comments: first, the CME in the 
language of CMW induces a charged dipole (of the vector density distribu- 
tion) that results from an initial nonzero axial charge density; second it has 
been recently pointed [3U] that an initial vector charge density via CMW will 
lead to a charged quadrupole distribution that may be observable in heavy 
ion collisions. For the rest of this contribution we will focus on the charged 
dipole signal for the CME phenomenon. 



1.2 Hunting for the CME in heavy ion collisions 

Now we turn to two key questions: can the Chiral Magnetic Effect occur in 
heavy ion collisions, and if so, what observables serve as unambiguous signals 
for the CME? 

The answer to the first question seems to be positive. Two elements are 
needed for the CME to occur: an external magnetic field and a locally nonzero 
axial charge density. The relativistically moving heavy ions, typically with 
large positive charges (e.g. +79e for Au), carry strong magnetic (and electric) 
fields with them. In the short moments before/during/after the impact of 
two ions in non-central collisions, there is a very strong magnetic field in the 
reaction zone [21] [13]. In fact, such a magnetic field is estimated to be of the 
order of « 10 18 Gauss (22] [23] (see also [24]), probably the strongest, 
albeit transient, magnetic field in the present Universe. The other required 
element, a locally non- vanishing axial charge density, can also be created in 
the reaction zone during the collision process through sphaleron transitions 
(see e.g. [16] for disucssions and references therein). As such, it appears at 
least during the very early stage of a heavy ion collision, there can be both 
strong magnetic field and nonzero axial charge density in the created hot 
matter. Therefore, the CME should take place, that is, an electric current will 
be generated either parallel or anti-parallel to the magnetic field B depending 
on the axial charge density is positive (due to sphaleron) or negative (due to 
anti-sphalcron). How large this current is, is of course another question, see 
e.g. [25]. 

The answer to the second question is much more difficult. Extracting the 
effects of the CME, which most likely occur at the very early stage of the 
collision, from the final observed hadrons, involves many uncertainties. First, 
it is quite unclear how long the magnetic field could remain strong: while the 
peak value is large, it decays very rapidly with time (if the only source of 
such field is from the protons in the ions) [26]. Second, if the CME current 
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is generated mostly at very early time, it is not clear to which extent this 
current could survive without significant modifications, since we know that 
the created quark-gluon plasma behaves like a strongly interacting fluid. Fur- 
thermore, even if this current survives, one has to find the right observable for 
its detection. At present, there is no satisfactory resolution on the first two 
issues. This will likely require comprehensive and quantitative model studies. 
In this review we will only focus on the third issue — the observables to 
be used for measuring the possible CME current and related "background" 
effects. 

In a simplistic view, one may consider the ultimate manifestation of the 
CME as a separation of charged hadrons along the direction of the initial 
magnetic field: more positive hadrons moving in one direction while more 
negative hadrons in the opposite direction. As a result, the momentum distri- 
bution of the final hadrons will have a charged dipolc moment. The direction 
of such a momentum space dipole is expected to be along the B field, parallel 
or anti-parallel, depending on the sign of the initial axial charge density in a 
given event. Since the initial axial charge may be positive and negative with 
equal probability, the event average of the momentum space dipole vanishes, 
< jy • B >= 0. This reflects the fact that parity is not broken globally by the 
strong interaction, so that any pseudo-scalar quantity, such as < jy • B >, 
will have to vanish. What one can hope for, however, is to measure the fluc- 
tuation or variance of this charge separation, i.e. < (jy • B) 2 >, which is a 
parity even quantity. As we will discuss later, the prize one has to pay is that 
other, conventional correlations, not related to the CME, may contribute to 
observables which are sensitive to the variance of charged dipole moment. 

Recently the STAR collaboration at the Brookhavcn's Rclativistic Heavy 
Ion Collider (RHIC) has reported [27] first measurements of a charge depen- 
dent correlation function in heavy ion collisions, which may by sensitive to 
the Chiral Magnetic Effect. The essential idea of the measurement, proposed 
by Voloshin [35], is based on two important features: first, in non-central 
heavy ion collisions, the direction of initial strong magnetic field is strongly 
correlated with the so-called reaction plane, which is spanned by the im- 
pact parameter and the beam direction. The B field is pointing (mostly) 
along the normal of reaction plane, albeit with random up/down orienta- 
tion; second, the CME-induced current, or the charged dipole in momentum 
space, implies particular charge-dependent correlation patterns. The same- 
sign charged hadrons will prefer moving together while the opposite-sign 
charged hadrons moving back-to-back along the B field direction, and thus 
perpendicular to the reaction plane, which is commonly referred to as the 
out-of-planc direction 0. While these measurements and their implications 
will be discussed in detail in Section [3] let us briefly summarize the present 



1 As a note of caution, the strong correlation between the B field direction and 
the participant-plane are considerably modified when the strong fluctuations in the 
initial condition are properly taken into account. As a result the two are rather weakly 
correlated in very central and very peripheral collisions | 23l 1241 129| . 
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status: the STAR (later PHENIX, and also ALICE) data show very interest- 
ing charge dependent azimuthal correlation patterns, and some features are 
in line with the CME predictions. Other aspects of the data, on the other 
hand, are very hard to understand within the framework of the CME. At 
present, therefore, the observation of the Chiral Magnetic Effect in heavy 
ion collisions, and the local parity violation in the aforementioned sense, has 
not been established experimentally, and additional measurements as well as 
further theoretical analysis arc required before definitive conclusions can be 
drawn. 

This review is organized as follows: in Section 2, we will present a gen- 
eral discussion on the charge-dependent correlation measurements in heavy 
ion collisions, with the emphasis on the CME related observables; in Section 
3, the presently available data from heavy ion collisions at a variety of col- 
lision energies will be examined and their interpretations will be critically 
evaluated; in Section 4, various possible "background" effects and their man- 
ifestation in various observables will be quantitatively analyzed; finally in 
Section 5 we summarize and conclude. 



2 The charge-dependent correlation measurements 

In this Section, we focus on various charge-dependent correlation measure- 
ments in heavy ion collisions and what can be learned from these observables. 
The emphasis will not be on the data themselves, which will be the subject of 
the next Section. Instead we will set up the conceptual framework for study- 
ing the azimuthal correlations, discuss possible complications in the design 
of the observables, and examine the connection between physical effects and 
the measurements. 



2.1 General considerations concerning azimuthal 
correlation measurements 

The basic experimental information about the (hadronic) final state of a 
heavy ion collision consists of the momenta and the identity - the electric 
charge, mass and possibly other quantum numbers - of all hadrons observed 
in the acceptance of a given experiment. Customarily, the three-momentum 
p is represented by the (longitudinal) rapidity, y, the transverse momentum 
p t as well as the the azimuthal angle <p. Events may further be grouped ac- 
cording to the charged particle multiplicity, which is a good measure of the 
centrality or impact parameter of a collision. From a given sample of events 
one can then extract the single particle distributions, d 3 N/dydp^d(j) either 



Charge-Dependent Correlations 



7 



for all charged hadrons or, more selectively, for identified pions, kaons, pro- 
tons, etc. In order to study possible correlations one analyses two-particle, 
three-particle and multi-particle distributions of various kinds. Most of the 
discussion in this review will focus on the dependence of various measure- 
ments on the azimuthal angle. The rapidity y and the transverse momentum 
Pt will either be in specific bins or integrated over. 

The analysis of azimuthal distributions has to deal with the fact that the 
azimuthal direction of each collision, characterized by either the direction of 
the angular momentum or the impact parameter, is randomly distributed 
in the laboratory frame. Therefore, a single particle azimuthal distribution, 
dN/dij) will always be uniform and, thus, rather meaningless. To learn some- 
thing about azimuthal distributions, one cither measures distributions of 
the difference of the azimuthal angles of two particles, dN/d(4>\ — ^2), or 
one determines the azimuthal orientation of a given event and studies dis- 
tributions with respect to this direction. Commonly the azimuthal direc- 
tion of the so-called reaction plane is used to characterize the orientation 
of an event. As already discussed in the Introduction, the reaction plane 
is spanned by the beam direction and the impact parameter of the colli- 
sion. Its orientation in the laboratory frame is given by the so-called reaction 
plane angle, &rp, which measures the direction of the impact parameter in 
the laboratory frame. Given the reaction plane angle, one then can study 
azimuthal angular distributions with respect to the reaction plane angle, 
f(<t> — &Rp) = dN/d((j> — \&rp). Clearly the determination of the reaction 
plane requires the measurement of other particles in addition to that used 
for the angular distribution (for a comprehensive review, see [30]). Therefore, 
the extraction of azimuthal distributions will require the measurement of 
two-particle (for the angular difference distribution dN/d(<f>i — 02 )) or even 
higher particle distributions. 

However, it is important to distinguish between the need to measure two- 
or many-particle distributions to study azimuthal distributions, and the pres- 
ence of true dynamical two- or many- particle correlations. To make this dis- 
tinction more transparent, it is useful to introduce an intrinsic frame or co- 
ordinate system where the x-direction is given by the direction of the impact 
parameter, which is typically referred to as the so-called "in-plane-direction" , 
and the y direction is defined by the angular momentum, or the so-called "out- 
of plane direction" . The relative angle of the x-axis of the intrinsic frame and 
that of the laboratory frame is then given by the reaction plane angle \Prp, as 
illustrated in FigHJ In theoretical considerations and model calculations the 
orientation of the reaction plane is assumed to be known, or in other words, 
these calculations take place in the intrinsic frame. Finally, the azimuthal 
angle ^ in the intrinsic frame is related to the laboratory angle 4> by 



$ = 4>- Vrp 



(3) 
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To continue, let us, as an example, consider a single particle distribution 
in the intrinsic frame 

h(<P) =fx{4>- Vrp) cx 1 + 2v 2 cos[2(<2>)] = 1 + 2v 2 cos[2(0 - & RP )} (4) 

which has an azimuthal asymmetry, characterized by the second Fourier com- 
ponent of strength v 2 . This kind of distribution, which will be relevant for 
the subsequent discussion, is important in the context of the observed az- 
imuthal asymmetries in heavy ion collisions, which are generally attributed 
to the hydrodynamics evolution of the system in non-central collisions. The 
parameter v 2 is commonly referred to as the elliptic flow coefficient. For a 
detailed discussion see |30j . The value for the elliptic flow parameter, v 2 , may 
be obtained by measuring the second moment of the angular distribution, 
(cos 2 (</> — ^rp)) . To this end we have to determine the reaction plane angle 
in each event, calculate the average moment in the intrinsic frame of each 
event and then average over events: 



cos [2 - <P RP )} 

events f ly W "\ 

E \m) E «»[2(fc-**p(0)]f (5) 



1 Agents , 1 



N, 



particle k—1 



In terms of the distribution function fi this can be expressed 

->„)» - l "^ M *-*ri m ^~ t " r)] <•> 

J air R p J - Wrp) 

Let us next consider the two-particle distribution 

/a(#i,*2) = /i(^i)/i(^2) + C(<Pi,<P 2 ) (7) 

where the first term is simply the product of the single particle distributions, 
and the second term, C (<?i,<?2) represents possible, true, two-particle cor- 
relations. Since the two-particle distribution depends on two angles, <Pi and 
^2, in general it will have terms which depend only on the difference of the 
angle ~ (<Pi — <P 2 ) = ((f>i — (f> 2 ), and which arc independent of the direction 
of the reaction plane. It will also have terms which depend on the sum of the 
angles, ~ (<Pi +<P 2 ) = {(f>i+4> 2 — 2\Prp) which are dependent on the reaction 
plane direction. This may be illustrated by inserting into Eq. ([7]) the single 
particle distribution, Eq. Q, and neglecting the correlation term, i.e., setting 
C(#i,# 2 ) = 0. In this case 



2 In reality the ability to express the actual measurement, as described in Eq. ([5|l . 
in terms of an average of moments of the intrinsic distribution over the reaction 
plane angle requires a detailed analysis of all non-flow effects and flow fluctuations, 
as discussed in detail in Ref. | 30 |. 
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/a(*i,*2) = /i(*i)/i(*a) 

- 2v\ cos[2(<?i - <P 2 )] + 2v\ cos[2(^>i + & 2 )] 

= 2v\ cos[2(</) 1 - fa)] + 2v\ cos[2(0 1 + fa - 2W RP )\ (8) 



The term ~ cos[2(0i — fa)] which depends on the difference of the angles can 
then be extracted by the measurement of the two-particle correlation 



The measurement of the term ~ cos[2(0! + fa — 2*Prp)] requires the deter- 
mination of the reaction plane, or at least a three-particle correlation mea- 
surement. For our example, Eq. (J5J), (cos[2(</>i — fa)]) ~ v 2 , and in fact this 
is one of the frequently used (and the simplest) methods for measuring the 
elliptic flow. However, this method suffers from the so-called "non-flow" [30] 
contributions, which are due to the correlation term we have neglected in our 
example. Our simple example also demonstrates a very important fact: sin- 
gle particle distributions, such as f\ do contribute to multi-particle azimuthal 
correlations. This will be essential for the subsequent discussion where one of 
the tasks will be to disentangle the effects from true correlations and contri- 
butions from the single particle distributions. 

The above discussion can be easily extended to three- (and more) particle 
densities with the same basic conclusions: 

• The n-particle density will have terms which do not depend on the reaction 
plane, and thus may be extracted by the measurement of appropriate n- 
particle correlations. It will also have reaction plane dependent terms, 




(9) 



Py(laboratory) 




Fig. 1 A schematic demonstration of the proposed simultaneous analysis of Qj and 
Q2 vectors in the same event. 
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which require the measurement of at least n+1 particle correlations or the 
determination of the reaction plane. 
• Unless not very carefully designed, multi-particle correlations will contain 
contributions from the single particle distribution. 

Finally, the measurement of angular correlations is of course not restricted to 
the second Fourier moment. Recently the harmonic moments, (cos[n(</>i — ^2)]), 
have been measured in order to study flow fluctuations [31]. These correla- 
tions may also be measured in a more selective way, such as correlations 
for particles with same or opposite electric charges (the charge-dependent 
correlations), correlations for particles with certain quantum numbers (e.g. 
baryon-strangeness |32j). or correlations for particles within or between cer- 
tain kinematic regions (e.g. the soft- hard correlations [33]), e ^ c - 



2.2 Measuring the charge separation through azimuthal 
correlations 

Let us turn to possible azimuthal correlation measurements as the signal for 
the Chiral Magnetic Effect. Specifically, as discussed in the Introduction, we 
have to find azimuthal correlations which are sensitive to a possible out-of- 
plane "charge separation" . 

We begin by defining what we mean by "charge separation effect" . Con- 
sider the distribution of final state hadrons in the transverse momentum space 
as schematically shown in the Fig. [T] If the "center" of the positive charges 
happens to be different from that of the negative charges, then there is a sep- 
aration between two types of charges which may be quantified by an "electric 
dipole moment" in the transverse momentum space. Such a separation may 
arise either simply from statistical fluctuations or may be due to specific 
dynamical effect such as the CME. We note that such a charge separation 
occurs already at the single-particle distribution level in the intrinsic frame. 
Let us, therefore, define a charge-dependent single-particle azimuthal distri- 
bution, which, besides a possible momentum-space electric dipole moment, 
also includes the presence of elliptic flow: 

f x ((j),q) ocl + 2v 2 cos[2(cj)-& R p)}+2qxdiC08((l>-&cs) (10) 

Here q and <f> represent the charge and the azimuthal angle of a particle, 
respectively. The parameters v 2 and d\ quantify the elliptic flow and the 
charge separation effect, while ^cs specifies the azimuthal orientation of the 
elcctric-dipole and &rp the direction of the reaction plane (see Fig. |T|). It is 
important to notice that an additional random variable x = ±1 is introduced. 
This accounts for the fact that in a given event we may have sphaleron or anti- 
sphaleron transitions resulting in charge separation parallel or anti-parallel 
to the magnetic field. Consequently the sampling over all events with a given 
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reaction plane angle, <?i?p, corresponds to averaging the intrinsic distribution 
f x over x, namely /i =< f x > x oc 1 + 2 v 2 cos(2</>— 2<?pp). Physically speaking 
this means that the charge separation (or electric dipole, being V-odd) flips 
sign randomly and averages to zero, thus causing the expectation value of 
any parity-odd operator to vanish. However, since < x 2 >= 1 the presence 
of an event-by-event electric dipole may be observable through its variance. 

For measurements related to heavy ion collisions one may reasonably as- 
sume particle charges to be \q\ = 1 which is the case for almost all charged 
particles, e.g., charged pions and kaons, protons, etc. We note, that the above 
distribution does not contain a directed flow term ~ cos(0 — ^Prp) for either 
type of charges, which is reasonable if the distribution is measured in a sym- 
metric rapidity bin. 

Finally one may also consider a pj-diffcrcntial formulation of the charge 
separation effect or charge separation effects associated with higher harmonics 
in the azimuthal angle <f>. We note that the charge separation term consid- 
ered in Eq. (|10[) is actually the lowest harmonic in a more general charge- 
dependent Fourier series expansion in terms of the azimuthal angle. Various 
higher harmonics may be present due to e.g. the occurrence of multiple topo- 
logical objects and their distributions over the entire transverse plane, the 
influence of transverse flow as well as the re-scattering of the CME current 
with medium. Here we concentrate the discussion on a possible measurement 
of the lowest harmonic that is most relevant to the CME current. 

Let us next discuss how the above defined charge-dependent intrinsic 
single-particle distribution contributes to the charge-dependent azimuthal 
correlations recently measured by the STAR collaboration in [37]. Note that 
here we are only considering the contribution from the charge separation term 
in Eq. (|10p . while there are certainly additional contributions from two- and 
multi-particle correlations which we will discuss later in Section 4. Specifically 
the STAR collaboration has measured the following two- and three-particle 
correlations [2~T] . 

(i) The two-particle correlation < cos(0i — <pj) > for same-charge pairs (++ 

/ ) and opposite-charge pairs (-1 — ). The contribution to this correlator 

due to the charge-dependent intrinsic single-particle distribution, Eq. (|10j) is: 

<*++/— =< cos(& - 4>j) >++/— = d\ (11) 
5 + - =< cos(4>i - <j>j) >+_ = -d\ (12) 

(ii) The three-particle correlation < cos(0i + <fij — 2<fi k ) > for same-charge 

pairs (i,j = + + / ) and opposite-charge pairs = H — ) with the third 

particle, denoted by index k, having any charge. The contribution to these 
correlators due to the distribution, Eq. (fT0|) . turns out to be 

< cos(4>i + <f>j - 2<j) k ) > ++/ — t k - an y = v 2 d\ cos(2 AWcs) (13) 

< cos(</) ?; + fa - 2<j) k ) >+_, k _ any = -v 2 d\ cos(2 A&cs) (14) 
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where "k-any" indicates that the charge of the 3-rd particle may assume any 
value/sign. We have also defined the relative angle of the charged dipole with 
respect to the reaction plane, A&cs = ^cs — ^rp- The purpose of correlating 
the charged pair with the third particle is to address the reaction plane 
dependence of the pair-distribution, as discussed in the previous Section, 
Sect. 12.11 Indeed, the STAR collaboration has demonstrated [27] that the 
above three particle correlator is dominated by the reaction plane dependent 
two-particle correlation function < cos(0; + fa — 2&rp) > and within errors 
they have found that 

< cos(fa + fa - 2cj) k ) >= v 2 < cos(fa + fa - 2& RP ) > (15) 



Based on the distribution Eq. (fTOj) we find the same relation between these 
correlation functions, since the reaction-plane dependent two-particle corre- 
lation is given by 

7++/— =< cos(fa + fa - 2& RP ) >++/— = d\ cos(2 A& C s) (16) 

for same-charge pairs, and 

7+ _ =< cos{fa + fa - 2V RP ) >+_ = -d\ cos(2 AVcs) (17) 

for opposite-charge pairs. 

To make contact with the predictions of the CME for the above correlation 
functions, let us assume for the moment that an accurate identification of 
the reaction plane could be achieved. In this case we may rotate all events 
such that ^S/rp = 0. Furthermore the CME predicts AWcs = t/2, and thus, 
for &pp = the charge separation term will take the form of ~ dsm(fa 
[TBI |2"T1 12"5] . If the only contribution to the above correlations would be due 
to the CME, a very specific pattern arises: 

7++/— =< cos( fa + fa - 2& RP ) >++/— = - d\ < 0, (18) 
=< cos (^ - 4>j) >++/-- = +d'l>0. (19) 

7+_ =< cos(</> ?; + fa - 2& RP ) >+_ = + d\ > 0, (20) 
=< cos(fa - fa) >+_ = - d\ < 0. (21) 

This pattern for the correlations 7 and <5, if seen in the data, would constitute 
a very strong evidence for occurrence of the CME in these collisions. However, 
as pointed out in [34] , and as we shall discuss in more detail in Section [3] the 
STAR measurements do not show the above pattern. For example, while 
in the above analysis for the same-charge pairs the correlators 7 and 5 are 
expected to be equal in magnitude but opposite in sign, i.e., 7++ = — <5++ the 



while 
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STAR data finds them approximately equal in magnitude but with the same 
(negative) sign. 

2.3 The vector analysis for measuring the charge 
separation 

When exploring an important phenomenon such a local parity violation, it 
is very useful to develop multiple observables which test its predictions, such 
as the Chiral Magnetic Effect. This is particularly the case in the present 
situation. The signals due to the CME are expected to be rather weak and 
the observables are not free from various backgrounds due to "conventional" 
physics, such as two-particle correlations. In addition, the interpretation of 
the STAR data is rather ambiguous. Therefore, it will be very helpful to have 
an alternative observable which is sensitive to a possible charge separation 
with specific azimuthal orientation. Currently there are a few proposals, for 
example the QI vector analysis [35], the charge multiplicity asymmetry cor- 
relations [36] , and the out-of-plane charge asymmetry distribution |37j . Here 
we focus on a detailed discussion of the Q\ vector analysis [35] . 

The QI vector analysis aims at a direct measurement of the intrinsic 
charge-dependent distribution in Eq. (fTU)) by identifying the charged dipole 
moment vector QI of the final-state hadron distribution in the transverse 
momentum space. The magnitude Q\ and azimuthal angle •f'f of this vector 
can be determined in a given event by the following: 

QI cos ^ = ^ cos fa 

i 

QfaniPf = £)giSiix& (22) 

i 

where the summation is over all charged particles in the event, with qi the 
electric charge and the azimuthal angle of each particle. This method is in 
close analogy to the Q\ and Q2 vector analysis used for directed and elliptic 
flow (see e.g. [30] )■ In the Q2 analysis one evaluates the charge independent 
quadrupole moment Q2 and its direction 1J2 in a similar fashion 

Q2 cos 2&2 = /J cos 2<j>i 

i 

Q 2 sin 2^2 = sin 2fa (23) 



The angle #2 is a measure for the reaction plane angle, Frp such that for a 
system with infinite many particles W2 — ^ ^rp- 
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Fig. 2 (Color online) (a) The Q1 and (b) \&3f\ distributions for the four different 
scenarios described in the text. 

Contrary to Q2, the charge dipole vector, Q1, incorporates the electric 
charge qi of the particles. The mathematical details regarding the observable 
Q\ and its relation to multi-particle correlations can be found in [35] . 

In each event, both angles \P£ and W% are determined from a finite number 
of final state hadrons (see Fig. [1]). While these angles correspond to their 
idealized expectations &cs and ^rp only in the limit of infinite multiplicity, 
their distribution and in particular the distribution of their difference, A& = 
&i — #2 should provide a good estimator for the magnitude of the charged 
dipole angle with respect to the reaction plane, A^cs = ^cs — ^rp- 

The combined Q\- and Q2- analysis will then provide distributions for the 
magnitude of the electric dipole, Q\, and its relative angle with respect to 
•^2, A&. This is demonstrated in Fig. [2] where we show the distributions for 
various scenarios calculated in a Monte Carlo simulation [35]. 

• The black triangles correspond to a "benchmark" scenario, where we have 
only elliptic flow but neither a charged dipole nor any true pair corre- 
lations. Therefore the resulting distributions for Q\ and \A^\ arise only 
from pure statistical fluctuations. 

• The red diamonds have been obtained by adding a physical dipole along 
the out-of-plane direction with a magnitude of d\ = 0.025, to the bench- 
mark scenario. 

• The green boxes are based on the case where back-to-back angular cor- 
relation for about 1% of the same-charge pairs but no dipole have been 
added to the benchmark scenario. 

• The blue stars result the case where same-side angular correlation for 
about 1% of the opposite-charge pairs but no dipole have been added to 
the benchmark scenario. 

As can be seen from the comparison in Fig. [5] and a more detailed dis- 
cussion in |35| , only the combined analysis of the distributions of angle and 
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magnitude, is able to distinguish between scenarios based on conventional 
two-particle correlations and those involving a true charged momentum space 
dipolc as predicted by the CME. As further discussed in [35] the final conclu- 
sion on the possible existence of an electric dipole will likely require a joint 
analysis of all three types of measurements, discussed in this Section: the Q\ 
distribution, the AW distribution, as well as the charge-dependent azimuthal 
correlations 7 and 5. 



3 Interpretation of the available data 

After having discussed the general aspects of charge dependent correlation 
functions in Section 2 we will now turn our attention to the actual measure- 
ments of such correlation function. Following the proposal by Voloshin [35] 
the STAR collaboration [27] presented the first measurement of the reaction- 
plane dependent charged-pair correlation function 

7 Q!(3 = (cosOc + - 2W RP )) (24) 

for pairs of particles with same, (a, /3) = (+, +), (— , — ), and opposite charge, 
(a, (3) = (+, — ). As already discussed in the previous Section, in order to ob- 
tain the correlator 7 aiJ g STAR measured three-particle correlation functions, 
and demonstrated rather convincingly that, within errors, they are related to 
the reaction plane dependent two-particle charged-pair correlation function 
by 

(cos(0 Q + - 2<f> k )) = v 2 (cos(0 Q + 4> p - 2<P RP )) = v 2 f a ,p (25) 

where v 2 denotes the measured elliptic flow parameter characterizing the 
elliptic azimuthal asymmetry. The results of the STAR measurement for 7 aiJ g 
are shown in the left panel of Fig. |3j 

Since the relation, Eq. Q25j) . has been established experimentally, we will 
concentrate our discussion on the charge dependent pair correlation function, 
la,p, Eq. (|24[) . Furthermore, we will choose a frame where the reaction plane 
angle is set to zero, Wpp = 0, so that 

7q,/3 = (cos(0 Q + (f>p)) . (26) 

In this frame the in-plane direction coincides with the x-axis and the out- 
of-plane direction points along the y-axis. Also the average direction of the 
magnetic field will be along the y-axis. 

Before we examine the STAR data more carefully let us recall what the 
prediction for the charge separation due to the Chiral Magnetic Effects are. As 
discussed in the previous section, the electric momentum space dipole induced 
by the CME will point (in an ideal situation) either parallel or anti-parallel to 
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Fig. 3 The data from the STAR collaboration for the reaction plane dependent 
correlation function (cos(0 a + <j>p — 2'Pb.p)) (left) and the reaction-plane independent 
correlation function {cos((f> a — 4>p)) (right) for like-sign and unlike-sign pairs. Also 
shown (lines) are results from various model calculations. The Figures are from |27| . 



the direction of the magnetic field, which in the frame where ^bp = points 
along the y-axis (neglecting fluctuations of the magnetic field [SSI) - Therefore, 
the charge separation due to the CME predicts to have pairs of same charge 
preferably moving together along the positive or negative y-direction. Pairs 
with opposite charge, on the other hand, are predicted to move away from 
each other along the y-axis. In terms of the azimuthal angles, (f> a , <f>p this 
means 



for same-charge pairs, and 



. 7T 7T , , 3tT 3tT. . _. 

W^T't) (27) 



,7r Sir. ,3tt it. . . 

(2, T )or( T ,-) (28) 



for opposite-charge pairs. Since 

,37T 37T. 

«*»(- + 3) = cos (y + Y )= (29) 

cos(^ + |) = 1 (30) 

the correlation function y a ,Pi Eq. (|2^|) . is expected to be negative for same- 
charge pairs and positive for opposite-charge pairs. While the STAR data, 
shown in Fig. [3l indeed show a negative value for same-charge pairs, the 
result for opposite-charge pairs is, at best, only mildly positive and, within 
errors, compatible with zero. Since opposite charged pairs are predicted to 
move away from each other, one may argue their (anti-) correlation should 
be weakened as these particles will have to traverse the entire fireball [T5] . 
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Therefore, at first sight, the STAR data may indeed show a first evidence 
for the charge separation pattern as predicted by the CME. However, the 
interpretation of the data is more difficult. 

The complication arises from the fact that the correlation function J a ,p 
does not unambiguously determine the angular correlation of the pair. To 
see this, consider a same-charge pair with angles (cj) a ,(j)p) = (0,7r). In this 
case the particles move away from each other in the in-plane direction. This 
is just the opposite of the correlation predicted by the CME, where the two 
particles are moving with each other in the out-of-plane direction. For both 
cases we get 

7T 7T 

cos(^> a ,(/)p) = cos(0 + tt) = cos(- + -) = -1. (31) 

Thus, the correlation function "f a ,p is not able to distinguish between same- 
side out-of-plane correlations and back-to-back in-plane correlations. How- 
ever, this ambiguity can easily be resolved by considering the reaction plane 
independent correlation function 

S a ,fs = (cos(0 Q - (j) P )) (32) 

which STAR has also measured, and we show their results in the right 
panel of Fig. [3] In the frame, where tpRp = 0, the two correlation functions 
may be decomposed in the in-plane ~ (cos(0 Q ) cos^^)) and out-of-plane 
~ (sin(0 a ) sm((frp)} components: 

7q,/3 = (cos(4> a + <j>p)) = (cos(0 Q ) cos(<j)p)) - (sm((/) a ) sin(0^)) , 

5 a ,p = (cos(0 a - 4>p)) = (cos(0 Q ) cos^)) + (sin(0 a ) sin^)) . (33) 

Qualitatively the STAR measurement in Au+Au collisions for both these 
correlation functions, ~f a ,p and d aj p for same-sign and opposite-sign pairs of 
charged particles, may be characterized as follows (see Fig. [3]): 

• For same-sign pairs: 

(cos(0 Q + <t>p)) same * (cos(0 Q - 4>p)) same < 0. (34) 
Using Eq. (j3"3")l this implies 

(sin(^ a )sin(0 j9 )) 8ome ~ 0, 

(cos(</> Q ) cos(0 /3 )) saTne < 0. (35) 

• For opposite-sign pairs we find that 

(cos(^> a + <f>p)) opposite 

(cos(^ a - ^)) > 0. (36) 
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Again, using Eq. (|33[) . this means 

{sm(4> a ) MM) opposite ~ ( C0S (M ^(M) opposite > °" ( 37 ) 

The decomposition of the actual data into the in-plane and out-of-plane 
components is shown in Fig. S) Obviously the correlations for same-charge 
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Fig. 4 Correlations in-plane (cos(if> a ) cos(0^)) and out-of-plane (sm(tf> a ) sin(if>p)) for 
same- and opposite-charge pairs in Au + Au collisions as seen in the STAR data. 



pairs are predominantly in-plane and back-to-back. This is exactly the op- 
posite of what has been predicted by the Chiral Magnetic Effect. This is 
illustrated in Fig. [5] were we have sketched the experimental situation for 
same-charge pairs based on the STAR data. For pairs with opposite charge, 
both in-plane and out-of-plane correlations have the same (positive) sign and 
magnitude. This implies that opposite-charged pairs move together equally 
likely in the in-plane and out-of-plane directions. This behavior can at least 
qualitatively be understood by resonance/cluster decays [36] or local charge 
conservation [38] . 

In addition to the data shown in Fig. |3] STAR has also analyzed the reac- 
tion plane dependent correlation function differentially as a function of 
the pair transverse momentum (sum and difference) and rapidity difference. 
Both these results are within qualitative expectations for a charge separation 
effect due to the CME [34] . Unfortunately, similar differential information is 
not available for the reaction plane independent correlation function, 5 a ^. 
Therefore a differential decomposition into in-plane and out-of-plane compo- 
nents, as we have done here, unfortunately is not possible at this time. Such 
information may help to further constrain possible background effects as well 
as predictions from the CME. 

Recently, the ALICE collaboration reported [39] the measurement of the 
same correlation functions for Pb+Pb collisions at a center of mass energy 
of yfs = 2.76 TeV, about ten times that of the STAR measurement. Just like 
STAR, ALICE determined the reaction plane dependent correlation function 
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Fig. 5 Schematic illustration of the actual STAR measurement (red) together with 
the predictions from the Chiral Magnetic Effect (black) for same-charge pairs. 

7 Qj ^ integrated over transverse momentum and rapidity as well as differ- 
entially. Within errors, the data for the integrated correlation function j a a 
agree with those of the STAR measurement, and the differential measurement 
show the same qualitative features. 

For the reaction plane independent correlation function, 5 a .p, on the other 
hand, the ALICE date differ from those by STAR. In particular, ALICE finds 
this correlation function to be positive for both opposite- and same-charge 
pairs. ALICE also provides the in-plane and out-of-plane pair correlations, 
(cos((j) a ) cos(cj)p)) and (sin(0 Q ) sinf^)) , respectively. Similar to the STAR 
measurement ALICE finds that for opposite-charge pairs the in- and out-of- 
plane correlations are nearly identical and positive. For the same-charge pairs, 
however, ALICE finds both in- and out-of-plane projections to be positive, 
with the out-of-plane correlation slightly larger than the in-plane projection. 
This finding would be in qualitative agreement with the expectations from 
the CME. Amusingly, early predictions [JUj for the collision energy depen- 
dence of the CME expected a smaller effect at the very high energies where 
ALICE has been measuring, largely due to the shorter duration of the mag- 
netic field. Of course the complex dynamics of heavy ion collisions and the 
various background contributions turn quantitative predictions for these cor- 
relation functions into a very difficult task, and a final resolution will require 
a systematic analysis of all available data at various energies. 

Given that the STAR data show an in-plane back-to-back correlation for 
same-sign pairs, one may wonder if there is still room for a charge separation 
effect due to the CME. This has been analyzed in [34] with the result that 
for the transverse momentum and rapidity integrated data, which the above 
analysis is based on, the backgrounds need to exactly cancel the CME induced 
charge separation. This may be just a coincidence, however. After all the 
data by ALICE show a different trend for the same-sign correlations. For 
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this situation to be clarified further differential data for the reaction plane 
independent correlation function, 5 a $ are required for both collision energies. 

Finally, as part of the RHIC beam energy scan program, STAR has mea- 
sured the reaction-plane dependent correlator ~f a ,p f° r various collision ener- 
gies [41] , see also [26] . They find that the difference for the correlator between 
same-sign and opposite-sign pairs decreases with decreasing beam energy. 
Such a behavior is expected from the CME. However, as we will discuss in 
the next section, all background terms scale with the elliptic flow parameter, 
V2, which is known to decrease with decreasing collision energy as well. 

To conclude this section, presently available experimental results concern- 
ing the CME are inconclusive. While the integrated STAR data disfavor the 
presence of the CME, the ALICE data allow for more positive conclusions. 
Clearly, progress requires data at lower energies as well as, and most im- 
portantly, differential measurements of both the reaction plane dependent 
and reaction plane independent correlation functions. In addition, given the 
rather unsettled state of affairs, measurements of other observables, such as 
the one proposed in the previous section, would be very welcome. 



4 Discussion of various background contributions 

As discussed in the previous sections the Chiral Magnetic Effect (CME) - if it 
exists - contributes to the reaction plane dependent two-particle correlator, 
first introduced in [55]. As in the previous section we denote the reaction 
plane dependent two-particle correlator by 



where <f>i and 4>2 are the azimuthal angles of two particles, and \Prp is the re- 
action plane angle. In the following we will distinguish between 7++/ , 7-1 

and 7 denoting respectively the correlator (|38[) for same-sign pairs, opposite- 
sign pairs, and the correlator without specifying the sign of measured parti- 
cles. 

As discussed in Section [5J a detailed measurement of 7 was performed both 
at RHIC by the STAR collaboration [37] and at the LHC by the ALICE 
collaboration [55] • However, as already discussed in the previous Section, the 
interpretation of experimental data is not straightforward since various effects 
can contribute to 7. 

The presence of elliptic flow allows for practically all "conventional" two- 
particle correlations to contribute to the reaction-plane dependent correlation 
function, 7. This can be easily seen from the decomposition of 7 into in-plane 
and out-of-plane projections, Eq. (|33p 



7 = (COS(01 + 02 - 2<f HP )) , 



(38) 



7 = (cos(^i) cos(0 2 )) - (sin(0i) sin(0 2 )) • 



(39) 
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It is quite obvious that even if the underlying correlation mechanism does 
not depend on the reaction plane it will contribute to 7 in the presence of 
the elliptic anisotropy v 2 . This can be seen in an extreme, though unrealis- 
tic, situation where all particles are produced exactly in-plane. In this case 
(sin(0i) sin(02)) = simply because there arc no particles in the out-of-plane 
direction and 7 = (cos(4>i) cos(</>2)). Obviously in this case, the presence of 
any two-particle angular correlation mechanism will result in a non-zero value 
of 7. 

In this Section we will focus exclusively on the contribution to 7 driven 
by the non- vanishing elliptic anisotropy v 2 . First we will derive the general 
expression which relates the elliptic anisotropy and the correlator 7 in the 
presence of arbitrary two-particle correlations. Next we will discuss a few 
explicit mechanisms that need to be understood quantitatively, before any 
conclusions about the existence of the CME can be made. In particular we 
will address corrections due to transverse momentum conservation (TMC) 
[42l l43l [44] and the local charge conservation [38] , both of which appear to 
contribute significantly to 7. In the last part of the paper we will discuss 
the possibility of removing, in the model independent way, the elliptic-flow- 
related background from 7. 



4.1 General relation 

In this part we derive the general relation between the elliptic anisotropy 
V2 and the two-particle correlator 7 in the presence of an arbitrary reaction 
plane independent two-particle correlations. 
By definition, the two-particle correlator 7 is 

_ / P2(0i,^2,zi,Z2,^i?.p)cos(</>i +<f>2- 2^ RP )d<f>id(j)2dxidx2 , 4Q , 
/ P2(<t > i,fa,xi,x 2 ,&Rp)d<f)id(f>2dx 1 dx2 

where, to simplify our notation, we denote: x = {jpt, rj) and dx = p t dptdr\. Here 
p t is the absolute value of transverse-momentum, while 77 is pscudorapidity 
(or rapidity) . p 2 is the two-particle distribution in the intrinsic frame with the 
reaction plane angle ^Ppp. It can be expressed in terms of the single-particle 
distributions, and the underlying correlation function C (see Section [5]) 

p 2 (<f>l, <f>2, Xl, X 2 , &Rp) = p((f>l,X 1 ,$ r R p)p(<p2,X2,1 , Rp)[l + C(<f> 1 ,<f>2,Xi,X 2 )]. 

(41) 

To simplify our calculation we assume the single-particle distribution to be 
P {^x^rp) = ^M[l + 2t- 2 (x)cos(20-2^p)], (42) 
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where p${x) and v 2 (x) depend solely oni= (puV)- We neglect higher mo- 
ments v n since their contribution to 7 turns out to be proportional to v n v m 
which is much smaller then the leading term ~ v 2 , see Ref. |45j . 

If v 2 (x) 7^ 0, the single particle distributions depend on the reaction plane. 
Therefore, the part of the two-particle density (|4i"j) involving the two-particle 
correlation function C depends on the reaction plane even if C depends only 
on 4>\ — 4>2- 

Here we want to concentrate on those correlations that depend only on 
A(j) = 4>i — 4>2 , namely the underlying correlation mechanism is insensitive to 
the reaction plane orientation. The correlation function may be expanded in 
a Fourier series 

a n (xi,x 2 ) cos (nA<f)) , (43) 

where a n (xi,x 2 ) does not depend on <pi and <p 2 . Substituting (|43[) and (|41[) 
into Eq. (|40|) . we obtain 

7=7^2 J Po(xi)po{x2)ai(xi,x 2 )[v2(xi) + v 2 (x 2 )]dxidx 2 , (44) 

where N = J po(x)dx 7 and we have assumed that a n << 1. 

Equation (|44[) explains why all correlation mechanisms with a non-zero 
ai [x\ , x 2 ) contribute to 7. For instance, it has been shown that transverse mo- 
mentum conservation (TMC) leads to a correlation function which depends 
on cos (A<f)) /Ntot where AT to t is the total number of produced particles. 
In this case ai (#1,2:2) oc l/AT tot . Let us also emphasize that all correlations 
that depend on the momentum difference between particles Ak = |ki — k 2 | 
also contribute to 7. In this case: 

C(Ak) =C{k\ + kl- 2kik 2 cos(Z\0)) , (45) 

which naturally leads to a non- vanishing a\ term. 

To summarize, Eq. (|44[) explains why transverse-momentum conservation 
[3H |43j [44] , local charge-conservation [38], resonance- (cluster-) decay [36] . 
and all other correlations with A<f> dependence contribute to 7. 



4-2 Transverse momentum conservation 

Very soon after publication of the experimental data by the STAR Collabora- 
tion it was realized that transverse momentum conservation convoluted with 
the non-zero elliptic anisotropy can lead to a substantial corrections for 7 
[4"T1 14UI |4"2] . This can be easily seen for the simplified situation where all par- 
ticles are measured (in the full phase-space), and where they all have exactly 
the same magnitude of transverse momentum |pi,t| = |pt|, i = 1, ...,-/V to t- 
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In the frame where "f^p = the correlator 7 may be written as 

. cosO, 

E^- 1 

or, alternatively, 



/S^te + frA (46) 



. = / (Si cos(^:)) 2 - (Ei sm(^)) 2 - E t cos(20 t ) \ ^ 

\ E^j 1 / ' 

where i and j are summed over all particles in the full phase-space. In the sim- 
plified scenario, where |p»,t| = |pt|, the conservation of transverse momentum 
implies 

J2 i cos(^) = ^. sin (<t>i) = 0. (48) 
Consequently we obtain 

/ E,cos(2^) \ -V2 
' \N tot (N tot -l)/ N tot ' V ' 

where iV to t is the total number of particles. Taking, for example, the central- 
ity class 40 — 50% we approximately have V2 ~ 0.1 and N to t ~ 1500 leading 
to 7 w —0.7 ■ 10~ 4 from TMC. This is roughly a factor 3 — 4 smaller than 
the experimental data for the same-charge pairs. This is only a simple esti- 
mation and a more realistic AMPT calculations [44] suggest that the TMC 
contribution is roughly factor 2 smaller than the STAR data. 

In a similar way we obtain for the reaction plane independent correlation 
function 

5 = (cosfoi J-. (50) 

JVtot 

In this case for N tot w 1500 we obtain <5 w —0.7 • 10~ 3 which is compa- 
rable or slightly larger in magnitude than same-charge data by the STAR 
collaboration. 

Similar results hold also in a more realistic situation, where only a small 
fraction of all particles is measured, and the magnitudes of transverse mo- 
menta are distributed according to the thermal distribution. This has been 
discussed in detail in [33], and we will only show the most important results. 

Using the central limit theorem and implying the global conservation of 
transverse momentum, the two-particle distribution function reads [46ll48ll43] 

/ \ I \ 1 \ ( 1 , 2 (Pl,x+P2,x) 2 {Pl,y +P2,y) 2 \ 

MPi,p 2 =s p pi )p (p 2 ) 1 + Tt \ -, 2 . - \' y 1 2 " J ■ 

y Ntot 2N tot (pi) F 2N tot {pl) F J 

(51) 
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where x and y denote the two components of transverse momentum. F de- 
notes that the appropriate average is calculated for all particles in the full 
phase-space. The single particle distribution, p(pi), is given by Eq. (|4"2"j) . 

Before we continue let us clarify one subtle point. Equation (|51[) is derived 
assuming that we first sample particles with a given V2 and next we conserve 
transverse momentum for all particles. In reality the opposite scenario should 
be considered. First we should sample partons/particles with a conserved 
transverse momentum, and after this the elliptic anisotropy V2 should be 
generated according to some dynamical model. Of course the second approach 
is much more challenging and renders analytical calculations difficult. At the 
end of this Section we will show that both procedures lead to comparable 
results for 7 and 8 and their transverse momentum distributions, however 
the rapidity distributions are quite different. We will come back to this point 
later. 

Using Eq. (|51l) and Eq. (|^0|) we can derive the following relations for 7 
and 8: 



and 



1 (p t ) 2 n 2v 2 ,n - V2.F - V2,F(v2.n) 2 
'N tot (j> 2 t ) F l-(h, F f 



1 (pt) n 1 + {v 2 m) 2 - 2v 2 ,f v 2 ,n 



Ntot {pDf 1 - (v 2 ,f) 2 

where we have introduced certain weighted moments of V2' 

V2 = -, — r , v 2 = r - 5 T • (54) 



(52) 
(53) 



(Pt) (P 



In the above equations F and Q denote averages that are calculated for 
all particles in the full phase-space, or for all actually measured particles 
in the restricted phase-space, respectively. Performing explicit calculations 
with reasonable assumptions about pt and r\ dependence of the single-particle 
distribution p(p t ,r]) and elliptic flow V2(pt,T)), we have found that for mid- 
central and peripheral collisions [43] 

7 ' N part « -0.005, 8 ■ N part « -0.05, (55) 

where N part is the number of participants, also referred to as wounded nu- 
cleons [4S] , 

To summarize, transverse momentum conservation results in a negative 
contributions to both 7 and 8, and they are of the same order of magnitude 
as the experimental measurement for like-sign pairs. More precisely they are 
a factor of 3 — 5 (very peripheral - mid-central) less in magnitude for 7, and 
a factor 1.5 — 4 (mid-central - very peripheral) larger for 8 than the STAR 
data for the same-sign correlator. While it is rather difficult to understand 
the data with only transverse momentum conservation it is interesting to 
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notice that the STAR experiment is sensitive enough to the effect of global 
transverse momentum conservation. 

Using Eqs. (|40I51[) we can easily calculate the dependence of 7 and S on the 
sum p + = (pi.t +p2,t)/2 and difference p_ = \p\j — P2,t|/2 of the transverse 
momenta of the pair. For the STAR data [27], 7 is growing roughly linearly 
with p+ and is approximately constant as a function of p_ . Interestingly a 
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Fig. 6 The two-particle azimuthal correlator (cos(</>i + <?!>2)} vs p+ = (pi,t + P2,t)/2 
(blue line) and p_ = \pi t ~ P2,t\ (red line) for N part = 100. The results are in 
qualitative and partly quantitative agreement with the STAR data for the same- 
charge correlator. Figure from Ref. |43 |. 



very similar behavior is obtained in the scenario with only global transverse 
momentum conservation. As seen in Fig. [6] the contribution of TMC to 7 is 
consistent with the data for p + > 1 GeV and underestimates the data for 
p+ < 1 GeV. As expected, for 5 very similar dependence on p + and p_ is 
obtained but rescalcd by a value of v-x- 

Recently, very similar results were obtained in the AMPT model calcula- 
tion [44] , where the transverse momentum is conserved on the event-by-event 
basis, and we will come back to this point later. 



4.2.1 Pseudorapidity dependence 

Since the contribution to 7 due to transverse momentum conservation is 
proportional to V2, one would naively expect that its (pseudo)rapidity de- 
pendence trace the rather mild (pseud)rapidity dependence of v-i. 

However, as shown in Ref. |42j , under quite reasonable assumptions we 
can obtain very similar rapidity dependence as in the experimental data. In 
the STAR measurement |27j the correlator 7 is maximum for \r\\ — 772 1 = 
and is approximately linearly decreasing to values consistent with zero at 
1 771 -772 1 « 2. 
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For simplicity of the argument let us assume that at the time t — all pro- 
duced partons/particles are distributed between two separate bins in rapidity 
with enforced global transverse momentum conservation. It means that if in 
the first bin the total transverse momentum equals Ki jt , in the second bin 
K.2.t = — Ki it in each event. Assuming further that all particles have the 
same magnitude of transverse momentum |p;,t| = |pt|, i = 1, iVtot, we 
obtain the following relation 

7 oc / cos (0») cos(<p k ) - Y sin(&) ^ sin(0 fe ) \ , (56) 

\iel k£2 iel fe£2 / 

and 

7 cx (K Xtl K x , 2 - K yA K yt2 ) = (Kl, - K^) , (57) 

where x and y denote the components of transverse momentum. It should be 
noted that here we calculate the two-particle correlator where one particle is 
taken from a bin number 1, and the second particle from a bin number 2. 

Now let us evaluate 7 at time t=0. In this case it is reasonable to assume 
that there is no elliptic anisotropy V2 and 7 = by definition. If we assume 
that bins are separated enough in rapidity so that there is no momentum 
exchange between two bins during the fireball evolution, i.e., the total trans- 
verse momentum Ki 4 is constant, then 7 = also in the final state, even 
if subsequently a non-zero vi is generated. Of course this simple argument 
demonstrates only that having the global TMC we can still obtain a non- 
trivial dependence of 7 as a function of r\\ — 772. In Rcf. this problem 
was studied in detail in the cascade model, where it was shown that the 
satisfactory description of the data can be obtained. 

We conclude that although the TMC probably cannot explain the data 
completely it gives rise to contributions which are of the same order of mag- 
nitude and which exhibit both the transverse momentum and rapidity de- 
pendence in qualitative and partly quantitative agreement with the STAR 
data. 



4.3 AMPT model 

It is desirable to study the correlators 7 and 5 also in an advanced Monte 
Carlo model, which allows for the conservation of transverse momentum on 
the event-by-event basis, and which generates reasonable values for the el- 
liptic anisotropy v%. Recently such calculation was performed in the AMPT 
model and the results are presented in Ref. [33] with the conclusion that 
the signal coming form AMPT is dominated by TMC. Indeed, the obtained 
results, summarized in Figs. [7] and [8j are in very good agreement with the 
previous discussion. As seen in Fig. [71 the calculated correlator 7 is in a good 
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Fig. 7 The two-particle azimuthal correlator 7 as a function of centrality in the 
AMPT model with different values of initial charge separation. 



qualitative agreement with the data, however it underestimates the STAR 
data by a factor of ~ 2. It was also shown that initializing the AMPT calcu- 
lation with the charge dipole leads to a better description of the data. 

In Fig. [5] the correlator 7 is plotted as a function of p + = (pi it +p 2 ,t)/2 
and Ar\ = rji — 772. Again, qualitatively the AMPT model reproduce the data 
but underestimate it by a factor of 2. 

As seen from the figures the AMPT model calculation (without initial 
dipole) is consistent with the global TMC and allows to understand the be- 
havior of the STAR data. It clearly demonstrates that all possible background 
effects must be studied very thoroughly before any conclusion about local 
parity violation can be reached. 



4-4 Local charge conservation 

Given the previous discussion, it is useful to construct a two-particle corre- 
lator which is insensitive to transverse momentum conservation and other 
charge independent correlations. The natural choice is the difference of 
opposite-sign and same-sign pair correlator see [55] : 



(58) 
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Fig. 8 The two-particle azimuthal correlator 7 as a function of p+ = (pi,t +P2,t)/2 
and Ar] = 771 —772 in the AMPT model with different initial values of charge separation. 

It is clear that only correlations that are charge sensitive will contribute to 7p. 
While the CME, if present, will contribute to 7p, global TMC, for example 
will not. Therefore the successful description of "fp with conventional physics 
would constitute a serious challenge for the interpretation of the data in terms 
of the CME. 

In Ref. [38j it was argued that 7^ can be fully understood assuming that 
charges are produced later in the collision (delayed hadronization). Indeed, 
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in the calculation of Ref. |38j the charges are produced in pairs in the same 
point in space-time. Due to the collective flow the initial correlation in space- 
time is translated into correlations in momentum space, and consequently it 
contributes to "fp. In this approach particles are emitted according to the 
blast-wave model with the additional requirement of local charge conserva- 
tion at freeze-out. Local charge balance is enforced within the finite range 
in rapidity a n and the azimuthal angle cr^. By comparing the model with 
experimental data on the balance function |50j . the values of and 0$ can 
be extracted which allows to make prediction for jp. 
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Fig. 9 Model with a delayed charge production and the local charge conservation vs 
the STAR data on the balance function in the relative azimuthal angle A<f> = <f>i — 4>2- 



The details of this calculation are presented in Ref. [38] . Here we summa- 
rize only the main results. It turns out that the model with delayed charge 
creation and local charge conservation can provide a successful description of 
the balance function both in the relative angle A(f> = 4>x — 4>2 and the relative 
pseudorapidity Arj, see Fig. |9]as an example^ 

Using the best values of parameters cr n and cr^ the contribution of local 
charge conservation to -fp can be calculated. The results are presented in Fig. 
M 

As seen in Fig. [TU]the agreement of the model with the STAR data is very 
good. It suggests that the two-particle charge sensitive correlations may be 
dominated by the local charge conservation. 



3 Recently a similar model was proposed to explain the fall-off of the same-side ridge 
in Arj [5J. 
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Fig. 10 Model with a delayed charge production and local charge conservation vs 
the STAR data on jp as a function of p+ and Ar\. 

4-5 Decomposition of flow-induced and 
flow-independent contributions 

From the analysis of the data in Section 3 as well as the discussion of various 
"background" effects, it appears rather plausible that the observed charge- 
dependent correlation patterns in 7 and S contain contributions from more 
than one source. In particular there are effects whose contributions to these 
correlations are flow- dependent, for example the transverse momentum con- 
servation (TMC) or the local charge conservation (LCC). On the other hand 
the CME, if present, is flow independent. Let us, therefore, attempt a decom- 
position of flow-induced and flow-independent contributions. 

We first consider correlation effects where the underlying correlation func- 
tion C is independent of the reaction plane orientation: 






x10" 3 
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C(0i,0a) ^ p{<t>u^Rp)p{<l>2^Rp)C{<t>i- <t>2), (59) 

where p is a single particle distribution. Note that the above is true for both 
TMC and LCC effects. A correlation effect of this type will contribute to the 
measured correlators as follows: 

la,p ~ v 2 F a> p , J Q ,^ ~ F Q!(3 , (60) 

with the factor F a ,p representing the strength of the effects, and (a, 0) is 

either + + / or H — . Both the TMC and the LCC follow this pattern 

albeit with opposite contributions. Thus F represents the total of all effects 
of this type. 

We should note, however, that the above relation, Eq. ([SO]) is a simplifica- 
tion of the exact relation between 7 and ^2(^,77), which is given in Eq. (pHf. 
Since the purpose of the present discussion is to gain some qualitative in- 
sight into the various contributions, we assume here that 7 is approximately 
proportional to the integrated v 2 ■ 

Next we consider possible contributions of the CME type. They would 
appear in the two-particle density in the following form 

P2(</>i,(h) ocsin^x - & RP ) sm(<f> 2 -&rp), (61) 

which explicitly involves the reaction plane. This term will contribute to the 
measured correlators as follows: 

7q,/3 ~ -Ha,p , <W ~ Hcc,p, (62) 

with the factor H representing the strength of the effects. It should be pointed 
out that besides the CME, there are possibly other effects that may also con- 
tribute to the correlators with the above pattern. One example is a possible 
dipole asymmetry from initial condition fluctuations that preferably aligns 
with the out-of-plane direction, see Ref. [52] for details. However, this effect 
will be charge independent, i.e., -£/++/__ = H + - > 0, whereas the CME 
predicts a charge dependence, = — £T+_ > 0. 

Combining the two types of contributions we arrive at the following de- 
composition for the reaction plane dependent and independent correlation 
functions, 

JaJ ~ V2 F a ,p — H a> p , 5 ai p ~ F at p + H a> p. (63) 

Given these relations, we may use the STAR data for 7q,^, S a p, and v 2 to 
extract the strength factors F ai p and H at p as a function of centrality. The 
result of such a decomposition is shown in Fig. [TT] 

Given this analysis we make the following observations: (a) Both com- 
ponents are charge dependent, i.e. there is significant difference between 
+ + / and H — ; (b) In both cases, however, the same-charge and opposite- 
charge signals are not symmetric with respect to zero. This may indicate that 
in each category there are likely more than one source of correlations; (c) 
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Fig. 11 The strength factors F a p (left) and H a ,p (right) extracted from the de- 
composition analysis (see text for details). The blue boxes and red diamonds are for 
+ + / and H — , respectively. 



There is a strong residual centrality dependence for both types component, 
although the dependence on centrality from v 2 has been removed. This may 
indicate that the correlations depend also on the multiplicity, which changes 
from central to peripheral collisions. 

Although, as already noted, the above analysis is qualitative, let us en- 
tertain a possible scenario, which would be consistent with the above ob- 
servations: The flow-induced signals may have two sources, the TMC with 

F l+/i- = F l MC < and the LCC with F ++f— = and F +- C > 
\FT_ C \ > 0. The flow-independent signals may be from two different sources, 
the CME with H ( ^ I / E _ > and H^ ME < and the dipole asymmetry from 
fluctuations (DAF) with H^fF = Hf AF > 0. Such a combination would 
indeed lead to correlations with magnitude and sign in qualitative agreement 
with the data. However, a quantitative analysis would have to be based on 
the exact decomposition based on Eq. (|44|) . Alternatively, one may attempt 
a separation of flow dependent and how independent contributions in exper- 
iment. How this could be achieved will be discussed in the following. 



4-6 Suppression of elliptic- flow-induced correlations 

In this chapter we will discuss the possibility of removing the elliptic-flow- 
induced background from the experimental data. 

As seen in Eq. (|44p all contributions due to correlations are proportional 
to the elliptic flow parameter, i>2- Therefore, it would be desirable to control 
or remove this contribution by a suitable measurement. There are essentially 
two ways to go about this. 
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First, as proposed in Ref. [S3], is to study collisions of deformed nuclei, 
such as U + U. By selecting very central, "face on face" collisions where the 
deformation of the nuclei is imprinted on the fireball, elliptic flow should be 
generated while at the same time the magnetic field will be very small. Should 
one observe correlations of the same magnitude and structure as ones already 
reported by STAR, this would identify their origin as being due to conven- 
tional two-particle correlations. The observation of considerably smaller cor- 
relations combined with a sizable V2 , on the other hand, would lend support 
for the existence of the CME. This approach, while challenging to analyze, is 
at present being attempted at RHIC, were first U + U collisions were made 
available. 

Alternatively, as proposed in Ref. [45] one may make use of the large cvent- 
by-event fluctuations of V2- By selecting events with different v 2 in a given 
centrality class we can control this background. In principle the measurement 
can be extrapolated to vi = (and consequently V2 {pt ,Tj) = 0) which will 
allow to extract correlations that only depend on the reaction plane orien- 
tation. Indeed, as presented in Fig. [12] even at b — 10 fm we expect a large 
fluctuation of initial eccentricity that will translate to large fluctuations of 
elliptic flow V2- 




Fig. 12 The distribution of initial eccentricity £2 calculated in the Glauber Monte- 
Carlo at the impact parameter b = 10 fm. 

Of course it is important to remove this background under the condi- 
tion that the contribution from the Chiral Magnetic Effect is approximately 
unchanged. In Fig. [13] wc demonstrate that indeed it is the case. Both the 
wounded nuclcons and spectators' contribution to the magnetic field weekly 
depends on €2- 

We believe this analysis should help to clarify the situation. Observation of 
non-zero 7++ or 7-1 at vanishing value of elliptic anisotropy V2 will suggest 
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e 2 

Fig. 13 The out-of-plane component of the magnetic field from wounded protons and 
spectator protons as a function of initial eccentricity £2 at a given impact parameter 
b = 10 fm. 



the existence of the correlation mechanism that is sensitive to the average 
direction of the magnetic field - possibly the Chiral Magnetic Effect. 

To summarize this Section, clearly there are many contributions based on 
conventional physics which contribute to the azimuthal correlations analyzed 
by the various experiments. In addition to those discussed in more detail in 
this Section other mechanisms, such as the decay of multi-particle clusters 
[36] . have been proposed. While it is more difficult to asses their contribution 
quantitatively, their influence cannot a priory be ignored. Therefore, it seems 
the best way forward is to separate the influence of elliptic flow and magnetic 
field experimentally as discussed in the last part of this Section. 



5 Summary and Conclusions 

In this review we have concentrated on the observational aspects of the search 
for phenomena related to local parity violation of the strong interaction. 
Specifically we have discussed various observables and their measurement for 
the charge separation, which is a predicted consequence of induced currents 
due to sphaleron and anti-sphaleron transitions in an external magnetic field. 
This phenomenon is often referred to as the Chiral Magnetic Effect (CME). 

We have discussed various properties and aspects of azimuthal angle cor- 
relations, and we have emphasized that, due to the elliptic flow observed 
in heavy ion collisions, virtually any two-particle correlations contribute to 
the azimuthal correlations. We have further discussed an alternative observ- 
able, which in our view may be better suited in discriminating between the 
backgrounds and the CME. 
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We have examined the presently available data on reaction plane depen- 
dent and independent correlation functions of same- and opposite-charged 
pairs. Our phenomcnological analysis of the data by the STAR collaboration 
showed that the measured correlations of same-charge pairs are predomi- 
nantly back-to-back, and in-plane. This is opposite to the predictions from 
the CME, where same-side out-of-plane correlations for pairs of the same 
charge are expected. The data by the ALICE collaboration taken at about 
ten times the STAR collision energy, on the other hand, show a correlation, 
albeit small, which is qualitatively consistent with the CME expectations. 

However, before any conclusion on the CME can be drawn, the contri- 
butions due to "conventional" correlations need to be accounted for. As we 
have discussed in some detail, both the conservation of transverse momen- 
tum as well as local charge conservation give rise to corrections which are of 
the same order as the experimental signal. These need to be understood and 
properly subtracted from the data in order to see if a signal consistent with 
the CME remains. Since there are conceivably many other two-particle cor- 
relations, which may enter due to the presence of elliptic flow, an important 
step towards answering the question about the existence of the CME is to 
experimentally disentangle the elliptic flow phenomenon from the creation of 
a strong magnetic field. This can be either done by colliding deformed nuclei 
or by carefully utilizing the fluctuations of the elliptic flow. 

In conclusion, the present experimental evidence for the existence of the 
CME is rather ambiguous. While progress on the assessment of the various 
background terms is to be expected, the sheer variety of possible correlations 
will likely limit a reliable quantitative determination of all the backgrounds. 
Therefore, the most important next step is the experimental separation of 
elliptic flow and magnetic field. In this context it is encouraging to note, that 
the first U + U collisions at RHIC have just been recorded. 

Finally let us close with a note of caution. Besides the CME there are 
other phenomena related to local non-vanishing topological charge fluctua- 
tions, such as the Chiral Magnetic Wave. One of the predictions in this case 
is the difference of elliptic flow between positively and negatively charged 
pions for collisions at lower energies [20] . However, again there may be other, 
more mundane effects which lead to similar phenomena, such as an increased 
stopping of baryon number and isospin at lower energies [54] . 
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